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We consider an unconventional superconductor in a low frequency harmonic magnetic field. In 
the Meissner regime at low temperatures a nonlinear magnetic response arises from quasiparticle 
excitations near minima in the energy gap. As a consequence various physical quantities acquire 
higher harmonics of the frequency of the applied ac field. We discuss how examination of the field 
and angular dependences of these harmonics allows determination of the structure of the energy 
gap. We show how to distinguish nodes from small finite minima ( "quasinodes"). Gaps with nodal 
lines, give rise to universal power law field dependences for the nonlinear magnetic moment and 
the nonlinear magnetic torque. They both have separable temporal and angular dependences. In 
contrast, with gap functions which only have quasinodes, these physical quantities do not display 
power laws in the applied field, and their temporal and angular dependences are no longer separable. 
We illustrate this via the example of the nonlinear magnetic moment for a, d + is gap. We discuss 
how to perform ac measurements so as to maximize the nonlinear signal, and how to investigate in 
detail the properties of the superconducting minima, thus determining the gap function symmetry. 
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I. INTRODUCTION 

There are strong-indications from numerous experimental results and theoretical calculations that the symmetry of 
the pairing stataJQ in various superconducting materials is not that of an isotropic s-wave. Unconventional pairing 
states have been assigned to diffe|ff;nt High Temperature Superconductors (HTSC's), Heavy Fermion33'El (HF) and 
some organic charge transfer saltajij (OS). Interest in determining the pairing state for these materials arises both 
from the efforts to obtain significant clues about the microscopic mechanisms responsible for superconductivity and 
to better understand their physical properties. 

For HTSC's it is widely accepted that most of the experimental results support a predominantly d-wave symmetry. 
There is however no consensus about the presence of admixtures of pairing states of different symmetry which would 
modify the position and value of the minima in the superconducting energy gap. These admixtures might cause the 
angle between the nodal lines to depart from 7r/2, or convert the nodes to very deep minima ("quasinodes"), or both. 
Part of the difficulties result from the surface character of many high quality pairing state probes. They measure 
information about the order parameter (OP) within a length scale of a few coherence lengths and are very susceptible 
to material imperfections near the surface |—|F|»rthermore it is not clear whether the pairing states are the same in 
the surface region as thixiughout the bulkJlSrlla There are also indications that the symmetry of the OP might be 
temperature dependent E3 113 The pairing state controversy for other suggested unconventional superconductors such 
as HF, SrRu204ll3't£l and some OS is even more ambiguous. 

In this paper we consider the low frequency magnetic response of a spin singlet unconventional superconductor in 
a time dependent magnetic field. We focus on the law-temperature regime in the Meissner state for OP's having lines 
of nodes (or quasinodes.) The nonlinear responsecErEa arising frorn-quasiparticle excitations near the minima in the 
energy gap generates in various physical quantities higher harmonicfS of the applied field frequency. Since the response 
extends over a length scale on the order of penetration depth, it constitutes a bulk probe of the superconducting OP. 
The use of the nonlinear response to a time-independent, field to perform gap spectroscopy (that is, to locate the 
nodes or quasinodes in the gap) was previously discussedtj. However, nonlinear effects are best detected through the 
use of ac techniques, since then the nonlinear response takes place at frequencies different from the input frequency 
at which the much larger linear response is found. These techniques significantly simplify the process of resolving the 
desired small nonlinear signal, which is a signature of the symmetry of the energy gap, from the large spurious but 
linear effects such as demagnetization factors, a — b plane penetration depth anisotropy, and trapped flux. 

In Section II we solve the nonlinear Maxwell-London equations in the low frequency Meissner regime. We generalize 
the perturbation method of Ref. ^ to include the temporal dependence. The method is illustrated on the example 
of OP with mixed d + s symmetry. We investigate the time and angular dependence of the nonlinear magnetic 
moment and the associated torque. The results for these quantities are easily extended to other forms of energy gaps 
with nodes since, as we shall see, in those cases one has separable temporal and angular dependences. The time 
dependence for gaps with lines of nodes is universal: Ha{t) \Ha{t)\ for the nonlinear magnetic moment and |-ffa(i)|^ 
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for the nonlinear magnetic torque, where Ha{t) is the appUed magnetic field. Both of these quantities have the same 
angular dependence. The nonlinear effects that we discuss here can be also viewed as field and angle dependent 
corrections to the superfluid density (penetration depth). We briefly discuss how our methods are suitable to extend 
studies of intermodulation and harmonic generationcj to low temperatures. 

In Section III we consider superconducting gaps without nodes but with quasinodes, as illustrated by a d + is OP 
with a small s component. We examine the nonlinear magnetic moment, which exhibits a more complicated temporal 
and angular dependence and a sharp enhancement of its maximum amplitude, compared to that occurring in the case 
of gaps with nodes. We show how to use these effects to experimentally distinguish nodes from small minima in the 
superconducting gap. In the final Section we present our conclusions and discuss possible extensions of this work. 



II. NONLINEAR MAGNETIC RESPONSE 



A. Maxwell-London Electrodynamics 

In the low frequency regime, i.e. in the quasi-static case|Hi~0 the relevant equations of the nonlinear Maxwell- 
London electrodynamics are formally the same as in the static case. Following the notation and results of the static 
or dc case,E3 we have 

VxVxv = ^j(v), (2.1) 

where the gauge invariant condensate flow field or superfluid "velocity" v is defined as: 

Vx , e 



V = 



2 c 



A, (2.2) 



with X the phase of the superconducting singlet OI^.^|t^e vector potential, and e the proton charge. The relation 
between j and v is generally nonlinear and given hy 



j{v) = -eNf d'.snis)^f[iwfv) + 2 dU{E{0 +^^f (2.3) 

JFS Jo 

where Nf is the total density of states at the Fermi level, n{s) the density of states at point s at the Fermi surface 
(FS), normalized to unity, v/(s) the s-dependent Fermi velocity, / the Fermi function, with E{^) = (^^ + |A(s)p)^/^, 



T the absolute temperature and A(s) the OP. The first term in Eq. (2.3) represents the supercurrent arising from 
the unperturbed condensate and the second is due to quasip article excitations. At T « 0, and for lines of nodes (or 
quasinodes), the second term of Eq. ( p.3| ) can be written a£3H 

j«p(v) = Y.^,pA^) « -2e^ A^/„ / ^^/«[(^/" ■ - " = 1'2, ... (2.4) 

n n "^^n 

where n labels the nodes (quasinodes) of the gap, f2„ denote regions where quasiparticle excitations are allowed, 
defined by |A((/)„)| -I- v/ • v < 0, is the azimuthal angle with ijespect to the closest node, n, and Vf « Vfn, its value 
at that node. Njn is the appropriate weighted density of states,E3 equ al to Nj for an isotropic FS. 

As stated in Ref. the relation between j and v given by Eq. ( ^.3[ ) can be used in the time dependent case 
provided tbat— the frequencies are smaller than the quasiparticle relaxation rate. For YBa2Cu3 07_5 (YBCO), this 
rate rangcsOLJ from 10^^ Hz to 10^^ Hz, dcpendmg on the temperature. Measurements of the low frequency nonlinear 
magnetic response are designed to be performedEj at frequencies under 100 Hz, well below this limit. 

For different OP's, we study the dependence of the nonlinear magnetic response on the direction of the field applied 
in the a — b plane. We considjet a slab sample, infinite in the a — b plane and of thickness d, much larger than 
the in- plane penetration depth,LZI in the c direction. This allows for an analytic solutiaa-and preserves the ang ular 
dependence of the nonlinear effects found for realistic finite three dimensional geometry.Ej'tj For a slab, Eq. (|]^) can 
be written as 



47re , , , , 

+ — j(v) =0, (2.5) 



and the boundary conditions are: 
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H(z,t) = HJi)|,^±a. 



(2.6) 



We take a particular time dependence of the applied field: 

Ha{t)^Hdc + HacCOSL0t, (2.7) 

which is suitable to experimentally study higher harmonics arising from the nonlinear response. The parity of quan- 
tities such as jqp under a sign change of Ha{t), as it occurs in an ac field, must be carefully taken into account. For 
an OP satisfying |A(— k/)| = |A(k/)|, where k/ is the Fermi wave vector, and for Ha chosen such that the nodal 
direction of quasiparticle excitations is £„, the key point is that a reversed field — Hq produces excitations along the 
opposite direction, — a;„. Therefore jqp is odd. This will help us to anticipate various modifications of the results for 
the static case. 



B. Gaps with nodes 

We consider an energy gap which can be accurately approximated near its nodesEf[E3 by 

|A(0„)| w |/xAe//0„|, n = l,2,.. (2.8) 
where /i is the slope of the dominant OP component and A^ff is the effective amplitude of the gap function. A 



particular case of Eq. (2.8) corresponds to a d ± s OP of the form Aj^±s{4>) = zLAg + A^i sin 20, where cj) is measured 
from the X axis, depicted in Fig. |l|. This form is often chosen to include the effects of orthorhombic distortion in 
the YBCO family of cuprates, with ^ = 2 and Ae// = (A^ - A^)i/^ For A^ ^ 0, the nodes of Eq. (U) are no 
longer separated by an angle of 7r/2. The nodal directions are along the unit vectors £„, n = 1,2,.. which form 
nonorthogonal axes. The nodes are shifted by an ang le ±v (see Fig. |l|), 

^^^sin-(^), (2.9) 

from the orthogonal axes, which we denote by X and Y (the ± sign corresponds to the d± s forms of the gap). 
Throughout this paper the X and Y axes will remain along the nodal directions of a pure d-wave. The direction of 
the applied field will be given by the angle between Ha(i = 0) and the +Y axis, depicted in Fig. |l|. It is convenient 
to introduce dimensionless versions of the superfluid velocity and magnetic field. For simplicity, we will perform our 
calculations below for an isotropic FS. Our results are largely independent of this assumption, as we will occasionally 
emphasize. It is however straightforward to extend these considerations and examine in detail the effects of FS and 
a — b plane penetration depth anisotropy as have been addressed in Ref. We define: 

u» = "'"^ , i^X,Y (2.10a) 

fiAeff 



Hdc,ac „ (^^eff , , 

— 77 — , -no — —; , n — rid 

Hq eXabVf 



(2.10b) 



where the in plane penetration depth, Aab, for a cylindrical FS is given by 2'itN fv ^ / c . In terms of these quantities 
and of the dimensionless coordinate Z = z/Aab, we obtain, from Eqs. (2.4), (2.f;) and (2.10): 



dzzUi - + sgn(/i(i))[eH('0)ux + e2i{^)uxUY + e3i('0)"i'] =0, i = X,Y. 



(2.11) 



where the factor sgn(/i(t)) arises from the odd parity of jqp with respect to h(t) = hdc + haccostot and the constants 
ei,2,3i('/') are defined in Appendix^. The appropriate boundary conditions for Eq. ( ^.11 ) are 



h{t) ^ , h{t) 

cosi/), Ozuy\z=z, ^ siu'tp, ux,y\z=o = 0, 



dzux\z=z 

where Zg = d/2Xab- The solution to Eq. (|2.1l|) and ( 2.12| ) can be sought in the form: 
u,{Z, t) = u^oiZ) + Y^KjiZ) cosjujt + u1j{Z) smjujt] j = 1, 2, 



(2.12) 



(2.13) 
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Using a perturbation methodHa, and taking into account the smallness of hac and h^c by neglecting cubic and higher 
terms in these parameters, we obtain the leading contribution to Ui{Z,t) (determined h>MLiQ, u'^-^ and The 
quantities of interest e.g., the magnetic mome nt or magn etic t orque , can be written downEa in terms of the surface 
values of the fields. By substituting Eq. ( 2.13 ) into ( ^.ll| ) and ( ^.12|) we get 



u.i{Zs,t) = h{t)[cosi>5ix + siiii>5.iY 
, h{t)\h{t)\ 



[eii cos V + cos -0 sin t/) + sin V'] = ^^Y, (2-14) 



where Sij is the Kronecker symbol. The angular dependence of the second term, nonlinear and nonanalytic in the 
field, is identical to that found in the case of constant applied field. Thus it remains to investigate the temporal 
dependence of various quantities which are also nonanalytic. The nonlinear transverse magnetic moment for gaps 
with nodes can be written (Appendix ^) as 

/, ,^ SXab Ha{t)\Ha{t)\ .3 3 

m± {^,t) = [eax sm'' -0 - eiy cos'' V 

D/iTT Hq 

+ cos?/'sinV'((eijf - e2y) cosi/' + {e2x - ear) sin-;/')] (2-15) 
^ iV„M^)iM/(v,) ^ N,^M^{i^,t) 

where S is the slab area and the factor Xab reflects that the nonlinear effects arise from currents present within a 
length scale on the order of penetration depth. We have introduced the normalization factor 

AT SXab{Hdc + Hacf 1 /o 1 «\ 

""-^ ^0 V! ^'-''^ 

and the normalized transverse moment AA^. The angular dependence of m±{ilj,t) or of M±{'4',t) is given by the 
function 

fit/j) = 3V3[e3x sin^ V - eiy cos^ V' + cos V'sin?/'((eix - e2Y + (e2X - ssy)) sinip]. (2-17) 

This result is independent of time and thus it is the same as in the static case. For a pure d-wave /(V") has a maximum 
value of unity, and so does the normalized moment as a function of angle. 

Experimentally, one can best detect the nonlinear effects by examining the dominant harmonics (time Fourier 
components) of m±{ijj,t), which are at 2uj and 3uj. Analysis of the field and angular dependence of these harmonics 
makes it possible to accurately determine the position of the nodes in the energy gap. One could obtain the same 
information from TO|| (-0, t), discussed in Appendix ^ but with greater difficulty because of the extremely large linear 
signal. The angular dependence of TO||(?/',i) is different from that of m±{^,t) but they both have identical temporal 
dependence. 

Since the measurements are performed in the Meissner state, one must have Ha{t) < where Hfi is the field 
of first fiux penetration, somewhat larger than Hd- Therefore, one wishes to optimize the experimental signal (i.e. 
the size of the 2lu or So; harmonics) , by an appropriate choice of the ac and dc components of the applied field at 
constant total maximum field. We can determine the optimal field mixture i.e., the ratio R = Hdc/Hac — hdc/hac (at 
fixed total field h = h^c + hac) which would produce the maximum signal for a normalized harmonic, Mj{tp) 

MM = r77^ r ^^^4r^ cos ju;td{u;t) 

cos ju;tdiu;t)f{^j) = M,/(i^), j - 0, 1, 2, ... (2.18) 



^(l + J,o)yo 

2 r h{t)\h{t)\ 



^(1 + '^io) 7o 

which factorizes into field and angular dependent parts, M j and fiijj). This factorization is also valid for an anisotropic 
FS. The field and angular dependence of m± for an anisotropic FS with nodal lines, as obtained in Ref. |2^, remains 
separable in a time varying applied field. As in the static case, FS anisotropy modifies the specific form of /{ip) but 
not the field dependence of the coefficient Mj , provided only that the characteristic field Hq is properly redefined. 
Therefore, the determination of the optimal R discussed below would remain the same. 

Since the angular dependence f{ip) is separable, determining the optimal field ratio R reduces to evaluating the 
coefficient Mj, the harmonic of the normalized angular amplitude of the nonlinear transverse magnetic moment. The 
normalizations chosen imply that in the static case, uj 0, Mq is unity. For j = 2, 3 we obtain from Eq. ( ^.18 ) 
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^ COS 375 h'^ 

M2 = -r^[hlj) + Ahdchac{cosp - {2hl^ + sm2p + sin4p] (2.19a) 

1 4 2 

M3 = ^[hl^ cos p - {-hl^ + -hlJ cos 3p + cos 5p ~ hdchac{2 sin 2p ~ sin 4p)] (2.19b) 

where p = sm~^ (hdc/hac)- Results for M2 and M3 at fixed maximum field as a function of the ratio R are given in 
Fig. H We show results for Hdc < Hac {R < I) only, since for larger values of R these harmonics are considerably 
smaller. The results indicate that the optimal applied field should be a pure ac field if the harmonic at 3w is measured, 
or Hac ~ Hdc for the 2uj harmonic. If one measures these harmonics for several values of R close to optimal, their 
dependence on R, as we shall see in Section III, can serve to experimentally distinguish nodes from small minima in 
the energy gap. The reduction of the amplitude for Mp-J = 2,3 compared to the static case (a factor of w 6) is more 
than compensated by the advantages of ac techniques.E3 

We can apply a similar analysis to the magnetic torque, t — mxtL. Torque magnetometry is an extremely sensitive 
experimental technique, which has already been used to investigate anisotropic magnetic properties of superconductors 
at higher fields.EjEd Here wc are interested in t^, which has an angular dependence identical to that of rn±{ip). The 
main harmonics arising from the nonlinear response are at Suj and 4w. Using a convention similar to that in Eq. 



( 2.18 ) we define 



2 r \h(f)\^ — 
TM) = - -j;r- cosju;tdiu;t)f{^j) = T,/(^), j = 3,4 (2.20) 

where we make use of the separability of the angular dependence. To determine the optimal R we need only consider 
Tj, the harmonic of normalized angular amplitude of the nonlinear torque. Elementary integration in Eq. ( 2.20| ) 
yields 

^3 = ^[f^p + 3hdchl,cosp-{h^a^ + 2hdchl,)cos3p (2.21a) 

3 3 sin4jD 
+ -hdchl^cos5p- {3hl^hac + ^hlJ{sm2p _il ) _ _^ sin 6p] 

1 /,3 POST) h'^ fi 3 

T4 = ^[^"^ - i^hlchac + -f )cos3p+{-hlKc + j^hlJ cos5p (2.21b) 

- ^ cos 7p - ^hdchl, sin 2p + [h^^ + ^hdchl,) sin Ap - ^^^^ sin 6p] 

The results for Tj as a function of R at fixed maximum field, seen in Fig. ^, show that for gaps with nodes it is 
optimal to measure the Auj time EC at pure ac applied field, and the 3a; component at i? w 0.35. As in our discussion 
of Mj , one can see that ES anisotropy does not alter Tj or the optimal value of R. 

C. Microstrip Resonator 

The nonlinear aspects we discuss here affect many other frequency dependent phenomena. We refer here briefly as 
an example to the nonlinear magnetic response of a microstrip resonator at low temperature. The particular geometry 
is that from Ref. ^ where the nonlinear dynamics was discussed for higher temperatures, when the relation between 
j and V is analytic. At low t emp erature the effect is quite different, as the nonlinearities are nonanalytic. For a 



cylindrical, isotropic FS, Eq. (2_^) can be expressed as: 

j = -epv + = -en^v, (2.22) 

where we use a notation distinguishing the electron density p from the superfluid density tensor fig , which also includes 
anisotropic, nonlinear effects due to quasiparticle excitations. The specific form of fig depends on the angular form of 
the energy gap and can be simply obtained from jqpi in Eq. (Al). The effects that we investigate can be viewed as 
arising either from a time dependent applied magnetic field or a driven time varying current. Following Rcf. |2^, we 
consider a pure d-wave gap with an applied field (or driven current) along nodal and antinodal directions. In both 
cases j II V. We have 
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A^(r = OJ = 0) _ n,(OJ) ]3l 

A2(0,j) -n.(0,0)-' ^'-'^^ 

where V' = corresponds (Fig. |l|) to j along the nodal and il) — 7r/4 to j along the antinodal direction. From Appendix 
A we obtain 6^=0 = 1/2 and 6^=^/4 = 1/2^/2. The critical current is jc = ns(0, 0)eArf/t;/ and in the absence of 
quasiparticle excitations at T = 0, j = 0, the superfluid density rig (0,0) equals p. 

We next examine the inductance per unit length, L, and show how the low temperature form of j(v) affects 
its nonlinear corrections. These considerations arc then readily applied to other quantitiGS-|Of interest such as the 
resistance per unit length or the quality factor Q. We use the expression for the inductancec3 

^ /,(g^ + A^(0,j)(47rj7c)^)dg 

with the integration being over a cross section of the microstrip having width w and thickness d. This expression can 
be rewritten as 

L = Lo + AL^ (2.25) 

where Lq is the inductance in the linear response, / = jgjdS and IcJgjcdS = wdjc- We make the approximatioi@ 
that the total correction AL due to the nonlinear response can be replaced by the kinetic part only, AL^i- Then, by 
using Eqs. ( |2.23D , ( p^ ), and ( p^ ) we obtain L ^ Lq + ALki \I\ / \Ic\ with 



47r6^A2(0,0)wd 




dS 
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AL^j = y-,^r^-Jsy (2.26) 



One can consider also the frequency dependent Q or the resistance, which can be decomposed into linear and nonlinear 
parts as in Eq. ( p.25| ). The power law behavior of the nonlinear corrections in these quantities i.e, linear in I/Ic, 
(instead of {iJJ.c)'^ at higher temperatures) will affect generation of higher harmonics. By using the transmission 
line equationsEj for a microstrip, one easily sees that these two types of nonlinearity give rise to second harmonic 
generation at low T eventually crossing over to third harmonic at higher temperatures. 

III. DETECTION OF NODELESS GAP FUNCTIONS 

In this section we investigate the quasi-static nonlinear magnetic response of anisotropic energy gaps with "quasin- 
odes" rather than with nodes. The nonlinear eiectromagnetism samples the .pairing state in the bulk. If the nodes 
are absent only in a layer of thickness much smaller than X^near the surfaceE3Ej, our results from Section ^ apply 
unmodified. We consider here a typical candidate, d + isJijCj energy gap with the d component being dominant .j-Tliej 
results presented are quite similar to those that can be obtained for other nodeless gaps such as d,j.2_y2 + idj^jjOcZl 
or an anisotropic s-wave. Generally, in the absence of nodes, the analysis of the low frequency nonlinear magnetic 
response is more complicated. 'T^Vtp are threshold effects that decisively affect the field dependence. 

We can approximate a d + zsOtj gap near its minima by 

|A(0„)|«(m^A202+a2)1/2^ n=l,2,.. (3.1) 
with As <C Ad and /i = 2 for the usual form of a d-wave OP. The small minima (A^) are located at the positions of 



the nodal points in a d-wave gap. The current due to the quasiparticle excitations is obtained from Eq. (2.4) 



j.p(v)^-:^E^/^/^" „2A2 -Q{v)(^r ■ x^f - Al) 71=1,2,... (3.2) 

Z LI L\j 



where the step function arises from the phase space for ahowed quasiparticle excitations (|A((/)n)| • v < ), given 
by 

<g < (3.3) 
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which is reduced compared to the d-wave gap. The step functions result in a threshold effect: it now takes a minimum 
field to create quasiparticles. However, as we shall later see, a small admixture of s-wave in a d + is gap can enhance 
the maximum value of the nonlinear effects at fields above threshold. From Eq. (3.3) it follows that the nonlinear 
effects are absent whenever Ha{t) is small enough so that \v\ < vt = ^s/vf- In dimcnsionlcss form the threshold field 
hx and the threshold velocity ut are 

hx = S, UT = —, S = -7^. (3.4) 



Substitution of Eq. (|3j) into (|2J) yields 

^zzu^~u, + sgn{h{t)){u^-u^^)e{\u,\-UT)^0, i = X,Y, (3.5) 

where Ui is given by Eq. ( 2.1C ) with Ae// = A^;. By modifying the methods used in this case for a static fieldil to 
include the time dependence we obtain 



Ui{Zs,t) = h{t)[cosipSix + sin^S^Y] + [^^(^) I (*^°^^ '^^^^ + ^^^^ i^^iv) - u'^sgn{h{t)) 

M\hit)\--, . J , ^^X,Y (3.6) 



3h(t) (cos ipSix + sin tpSiy ) (cos ipSix + sin ipSiv ) ' 

for the solution at the surface of a slab, Z = Zg {Zg 3> 1), which suffices to express the quantities m± and r^. 
Nonlinear effects are present only for h{t) > hTi{iJ,S) > hx, where hxx = S/ cosip and Hty = S/ sin^p are the angle 
dependent threshold fields required to excite jets of quasiparticles centered along the X and Y axes respectively. The 
transverse magnetic moment is 

m^{^,t) = ^ ^a{t) \Ha{t)\ cog^gjnV'[cosV'e(|/t(i)| " hrx) ~ sm^pe{\h{t)\ - hrv)] 

+ SS^Ho sgn(ffa(i))[cos Ve(|/i(t)| - hTv) - sin V'e(|/i(i)| - hrx)] (3.7) 
- 2<53-^[cot7^e(|/i(i)| - hry) - tanV'e(|/i(t)| - hrx)]], 

where we see that, as stated above, the temporal and angular dependence of m±{t,'il)) are no longer separable. To 
investigate the behavior of this expression at different fields and admixtures of s-wave OP it is useful to employ the 
scaling relation: 

m±{Kh{t), kS) = K^m±{h{t), S) (3.8) 

valid for any angle ip, at fixed field ratio, R. The quantity n is an arbitrary scaling factor. 

Again, we investigate the dominant harmonics of m^ (ip, t). Because of the non-separability of temporal and angular 
dependences, one cannot introduce M j as in Eq. (|1|). One must now carefully consider both the harmonics 

MjW = - -"7' ^ cosju;td{u;t), j = 2,3 (3.9) 

and the relevant angular Fourier components (EC's) in terms of which the angular dependence can be analyzed. The 
most important angular EC is that which reflects the main angular periodicity, 7r/2 of the energy gap 

Mf = -I Mj{tP) sin AiPdip, j = 2,3 (3.10) 



We will focus in the rest of this section on the dominant response of the nonlinear magnetic moment at 2uj and 30;. 

In Eig. ||we show our results for the angular dependence of the 3a; harmonic, M3 (?/;). The details of its calculation 
are given in Appendix We display the range ^ g [0, 7r/4]. We have taken R = (guided by Eig. |^) and a 
dimensionless field amplitude h — hac — 0.05. This value, using typical material parameters for YBCO, corresponds 
to Hac ~ Hfi. The solid line represents the d-wave result, normalized such that the maximum of Mj_{'ip) for the 
static case is unity. This normalization is employed throughout this section. The broken lines are labeled by the 
corresponding ratios S. 
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An obvious feature seen in Fig. ^ is the enhancement of the maximum value of the signal that may arise from small 
admixtures of an s-wave OP. At fixed field, the maximum enhancement of M^^ip) is much more pronounced than the 
moderate one for ■m±{ip) found in the static case (Fig. 7 in Ref. |2^). A similar enhancement is also present i n M -ijxj}). 
This result is rather unexpected since such an admixture would reduce the available phase space (recall Eq. (3^)) for 
quasiparticle excitations. This enhanced signal is of non-negligible significance in planning experiments. We therefore 
pause here to explain its physical origin. 

The nonlinear transverse magnetic moment arises from contributions due to components of the quasiparticle current 
jgp,n(v) at different nodes (quasinodes). These contributions partly cancel each other as each node, in effect, tries to 
twist the magnetic moment away from itself. Thus, if a magnetic field is applied along an antinode, the cancelation 
of such components is complete and mj_ = 0. To explain the peculiar enhancement in M^^ip), we need to investigate 
how the introduction of increases the asymmetry of such contributions from different jets. The asymmetry, which 
is reflected in the angular asymmetry of the curves in Fig. p] for (5 > 0, is brought about by one of the nodes being 
above its angular dependent threshold field (see below Eq. ( |3.6| )) while the other is still below. It was seen in Ref. 
p5| that this was responsible for the small enhancement seen in the dc case. Now we will see quantitatively how the 
larger enhancement in the ac case arises from the same source. 

We denote the contributions to A^j^ = m±/Nm of the individual jets along the X and Y axes at a fixed angle ipo 
by M±{ipo,t)x, M±{ipo,t)Y: 

M^{i'o,t) = M±(iJo,t)x + M±{i^o,t)Y. (3.11) 

In Fig. ^, computed for the same fields as in Fig. ^, we show A4±{ijj,t)x , ■M±{ip,t)Y and their sum, A4j_{il;,t), as a 
function of time ujt at a fixed angle ifj = %Ijq = 25" in the region where the maximum of M3{ip; 6 — 0.02) occurs. We see 
in Fig. U that the enhancement in this case is large, a factor of about two over the d-wave signal. The quantities in the 
figure are shown for both 6 = and 0.02. The curves are labeled by two indices, the first a letter denoting the jet, or 
sum of jets, plotted, and the second being the value of S. The range cut £ [0, 90°] shown in the figure could be extended 
over the whole period by symmetry. For the chosen value of V'o, the threshold fields satisfy hrxii^o) < ^Ty(V'o) and 
consequently the magnitude of the jet along the X axis is greater than that of the one along the Y axis. With the 
admixture of an s-wave OP, there is a region, tot ^ 20", where the magnitude of Ai±('4'o,t)x is reduced less than 
that of A4±{ipQ,t)Y- Therefore the cancelation of the contributions from the two different jets takes place to a lesser 
extent, resulting in an overall increase of A4±{^lJo,t). However this accounts for only part of the enhancement seen 
in Fig. ^, basically the part found in the dc case, which corresponds to that at cji = in the Figure. The remaining 
part of the increase in M^{ ip) is explained by the convolution of the time dependence of Mj_{'>Jj,t) and cosSujt in the 
definition of MJijj), Eq. (3.9). When the curves representing ■m±{ip,t) for 6 = 0, 0.0 2 an d cos3ujt, (which is also 
shown in Fig. ^ are multiplied the results obtained, which are the integrands in Eq. (3.9), are as shown in Fig. 



The solid curve represents the pure d-wave result and the broken curve that for 6 = 0.02. The ratio of the two net 
areas between the curves shown and the horizontal axis is about a factor of two, i.e. the factor between M3{ipo) for 
6 = and for 6 = 0.02 in Fig. ^. Thus, the enhancement is explained. 

It is straightforward to perform a very similar analysis for M2{ip) using its explicit expression from Appendix |^. 
From parity considerations, similar to those made for gaps with nodes, M2{ip) vanishes identically as a function of 6 
and ijj for h^c = 0. In this case, one chooses R « 0.5 (see Fig. ^) and again finds an enhanced signal, as in Fig. |^. 
The size of the enhancement is largely due to the convolution with cos 2ujt. 

The non-separable temporal and angular dependences of m^{ip,t) require additional care in trying to determine 
the optimal mixture of dc and ac fields in Ha{t) which maximizes the harmonics and angle EC's of the signal. From 
the expressions for A'hj ip) in Eq. ( plj ) we have obtained results for its 7r/2 angular FC. We denote this quantity by 
Af|, as defined in Eq. ( 0.1C| ). We show typical results in Fig. ^ There we plot M| as a function oi R a,t h = 0.05, 
for several values of 6. The solid line represents the pure d-wave and the broken lines are labeled by their value of 6. 
We see that the optimal field ratio depends on the particular s-wave admixture. Through the scaling relation in Eq. 
(3.8) these results can be simply extended to other values of h and 6. In Fig. ||, we plot M| (i.e. the tt/2 angular FC 



of the signal at 3a;) as a function of R at maximum dimensionless field h = 0.05 for the same values of 6. In this case 
R = produces the maximal signal, for all 6. 

It is instructive, in order to learn how to distinguish nodes from quasinodes, to contrast Fig. & which displays 
Mj, j = 2, 3, the angular amplitude of the 2uj and 3uj harmonics for a gap with nodes, with Figs. and |[ Because 
of the separable temporal and angular dependence of m^{ip,t) when nodes are present, the solid lines (pure d-wave 
case) in Fig. ^, |^ are simply proportional (with a normalization factoio ~ 1) to the curves M2 or M3 respectively. 
Experimentally, a sensitive test would be to measure , Mg for several values of i?, and compare the R dependence 
with the corresponding solid lines in Fig. The deviation from these lines would indicate absence of nodes in the 

gap. A further way of analyzing such data is suggested by Figs. ^ and |l^ where we sho w and M| respectively as 



functions of 6, for several values of R (additional results are easily obtained using Eq. (3.1C) and the expressions in 
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Appendix C). The curves plotted are labeled by their values of R and are given at a fixed maximum field h = 0.05. 
To compare experimental data with an assumed d + is gap, and to determine the s-wave admixture, one would look 
for a value of 6 which agrees with a particular dependence of the measured M| or M| at different R. In these figures 
we notice again the implications of the non-monotonic behavior of m± (or M±) with 5. Just as we found above for 
the harmonics, there is also a substantial enhancement for the angular FC at given R as compared to the S — value. 
It has the same origin, as shown in Fig. |^, ^, predominantly due to the convolution of A/j(?/', t) and cos{jLLjt), j = 2, 3. 
In Figs. ^ and |l^ we see that maximum such enhancement (at different value s of R) is located typically in the range 
0.015 ^ S < 0.025 at /i = 0.05. More generally, using the scaling relation Eq. ( |3.8| ), for an arbitrary maximum field h 
(within the Meissner regime) the maximum signal for Mj, j = 2,3 is approximately in the range h/3 ^ S ^ h/2. 

Thus, there are several distinct features for which the results for a nodeless gap with small minima differ considerably 
from those for gaps with true nodes. These results suggest how to perform measurements and analyze the nonlinear 
quantities in the low frequency, low temperature regime. They can be supplemented with the results from the static 
case. It would also be useful to study the field dependence oim± and its Fourier components at fixed R while changing 
the maximum field {Hdc + Hac)- 

Using the methods that we have presented, similar results for the behavior of the nonlinear magnetic torque can 
be obtained. Analogously, results for other gapped states such as d^2_y2 + idxy can be derived without difficulty. 

IV. CONCLUSIONS 

We have examined properties of the nonlinear magnetic response of an unconventional superconductor at low 
frequencies. There are two key points to our conclusions. First, these frequency dependent nonlinear phenomena 
provide a particularly powerful tool to distinguish nodes from quasinodes. Second, the experimental sensitivity that 
allows one to perform "node spectroscopy" can be considerably enhanced over the corresponding dc result. 

For gap functions with nodal lines the nonlinear moment (m^, m||) and torque (r^) have separable field (and hence, 
time) and angular dependences. These quantities are respectively quadratic and cubic (up to a sign determined by 
parity) in the instantaneous applied field. The nonlinear magnetic response for nodeless gaps is more complicated: 
the time (or field) and angular dependences are not separable and the field dependence is not a simple power law. 
We have discussed how these features can be studied by measuring higher harmonics (referred to the frequency of the 
applied ac field) of various physical quantities. For an applied field Ha{t) — Hj^^ + HacCOSLut, at a fixed maximum 
Hdc + Hac we have computed the optimal ratio R = Hdc/ Hac at which the signal (for a certain higher harmonic) 
is largest. For gaps with lines of nodes it is most favorable to measure the 2uj harmonic of m± a,t R ~ 0.5 or the 
3uj harmonic at i? = 0, while for the best alternatives are the 3w harmonic at i? w 0.35 and the 4a; harmonic 
component at i? = 0. These results are independent of the specific angular dependence of the gaps with nodal lines 
and are also valid for an anisotropic FS. In the absence of nodes, as we illustrate through the example of a d + zs gap, 
selecting the optimal value of R to maximize higher harmonics is less simple but more rewarding. The result depends 
on the particular admixture of s-wave component to the OP and hence provides for the quantitative determination 
of this admixture. The results we shown are sensitive to the angular position of nodal (quasinodal) lines, allowing for 
detailed characterization of the minima in the energy gap. 

The analysis presented here can be generalized to include point nodes, and to triplet pairing. The field dependence 
for gaps with point nodes is different than that found for gaps with lines of nodes. For example, m± is cubic rather 
than quadratic in the instantaneous value of the applied field. Using these considerations and the temperature range 
attainable by a dilution refrigerator it should be feasible to experimentally study other candidates for unconventional 
superconductivity, such as Heavy Fermions, Organic Superconductors and SrRu204. Analysis of the material pa- 
rameters (to estimate Hq and -ff/i) shows that in many cases the expected signal should be at least as large as for 
HTSC's. Experiments-investigating magnetic response in these materials, at sub-Kelvin temperatures, are currently 
under considerations .Ea 
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APPENDIX A: CURRENTS 

Integration of Eq. (2^) without a — b plane FS anisotropy yieldJii 



jqpi ^ o^Nf / [eui^)vjc + e2^WvxVY + e3i(V')4]: i = X,Y, (Al) 

2 nAeff 

where the coefficients ei^2.3i('0) S'^e defined for various ranges of tjj. Depending on the particular orientation of Hq, 
as given by the angle ip, these coefficients describe contributions from different jets. The results are quadratic in the 
superfluid velocity. For a d ± s OP we obtain: 

eix = ear = - , e2x = Sesx = 2eiy = e2Y = -2cs(c - s), ^ € [i^, f ~ ^1 (A2a) 

eix = -e^Y = -{(? + s^), e2x = -2e3x = 2eiy = -e2y = 2cs{c + s), G - z^, tt + i/]. (A2b) 
where c = cos i', and s = sinz^. The pure ci- wave limit corresponds to c = 1 and s = 0. 

APPENDIX B: MAGNETIC MOMENT FOR A GAP WITH LINES OF NODES 

The magnetic moment can be calculated using only the surface values of the fields HSS For the geometry con- 
sidered here, the magnetic moment, 



^y"drrxj(v), (Bl) 



2c 

can be expressed as: 

Sd Ha X V , , , /T-,r^\ 

where x and y are orthogonal axes fixed in space, S is. the slab area, and we have used that v is odd in z. For a slab, 
{d 3> Aafc), which is the case of experimental interest,Ej the nonlinear magnetic moment can be obtained from Ui{Z, t), 



given in Eq. (2.14). For a gap with node lines it is oc h(t) \ h(t) \. Th e transverse component, — mx cosV-'+m-y sin'0, 
perpendicular to the applied field, is obtained from Eq. ( |2.15 ) and the longitudinal component of the nonlinear 
magnetic moment, m|| — ~mx sinip + my cosip, is given by 

_ SXab Ha{t)\Ha{t)\ ■ 3 ; , 3 ; fT3o\ 

FT eix sm'' V + e3y CDS'" V (B3) 

+ cost/i sin?/;((e2x + siy) cosip + (e3x + e2y) sinV')]. 

For a pure d-wave OP and -0 S [0,7r/2], the angular dependence of m± (its normalization as given in Eq. ( ^.17[ )) 
is f{ip) — 3-^/3 cos V' sin -0 [cos V' — sin-0] and that for m|| it is f\\{ip) — 3-\/3[cos^ -0 + sin'^i/']. While the angular 
modulation for m|| is somewhat smaller than that of mi. the maximum of r7i|| is enhanced by a factor of 3\/3. This 
enhancement combined with ac measuring techniquesES could facilitate examining the higher harmonics due to the 
nonlinear response. 



APPENDIX C: FOURIER COMPONENTS FOR A NODELESS GAP 



W e give here the analytic expression for the harmonics Mj{tp), valid for a d + is gap and hdc < hac- Integration of 
Eq. ( 3.10 ) and substitution of m± from Eq. (3/?) yields 



2 3\/3 

MjOj) = — [cos0sinV'[cos0(pijei +^3^63) -sinV'(p2j62 +^4^64] 



TT /l2 

+ 3(5^[- sin -0(91^61 + g3j©3) + cosij{q2jQ2 + 94j©4] 
- 26^[- tan0(rijei +r3je3) + cot0(r2je2 + r4je4]]. 



(CI) 



J = 2,3 
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where we have introduced the abbreviations 

61,2 EE e[a + - /iTX,y (V', S)], 63,4 = e[b-a- hrx.Yii^, S)], 



(C2) 



and a, b are replacing h^c and hac respectively. Since we work at a fixed maximum field {hdc + hac — cons t.), and 
R = hdc/hac, R e [0, 1), we set a = {hdc + hac)R/{l + R),b= [hdc + hac) /{I + R)- The coefficients in Eq. for 
the 2a; component are 



Pj2 = — Wi + 2a6sinwj + [a? H ) sin2u'i H — sinSw^ H sin4u'i, i = 1, 2 

2 2 3 8 



(C3a) 



52 ^2 2 52 

Pfc2 = -rri^k - tt) + 2a6sinwfe + (a^ + — ) sin2wfe + - sinSw^ + — sin4wfe, fc = 3, 4 

Z Zoo 



(C3b) 



qi2 = sin 2wi 



i = l...,4 



(C3c) 



n2 



■ sm w,- 



b ' 62 
Those for the 3w component are 

62 



4a 2(2a2/52 _ 1) 

+ In 



(62 -a2)i/2 



(6 -a) tan^l^ + (62 -a2)i/2 



(6 -a) tan^l^ - (62 -a2)i/2 



Pi3 = "TT si'^ ""^j + si'^ H sm Swi + — sm Awi 



gi3 = -sin3wi, i = 1, ..,4 



= TrK^fl^ ^ b^)wi — Aabsmwi + 6^ sin2wi] 
6'' 
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■ sin bwj 



2a(-4a2 + 362) 
63(62 - a2)i/2 



In 



(6 -a) tan2|^ + (62-a2)i/2 



(6 -a) tan^ - (62 -a2)i/2 



1,..,4 



1,..,4 



1,2 



(C3d) 

(C4a) 
(C4b) 
(C4c) 



r-fcs = ^[(4a^ - 6^)(7r - Wk) + 4a6sinii;fc - 6^sin2wfe] 

(6 -a) tan^f + (62-a2)i/2 



63 

2a(-4a2 + 362) 



63(62 -a2)i/2 



In 



(6 -a) tan^f -(62-a2)i/2 



k = 3,4 



(C4d) 



where the ?«i,2,3,4 denote possible values of wt for which |/i(i)| = hTx.y- The integral in Eq. ( |3.9D can be replaced 
by integrals over intervals € [0, Wi], i = 1,2 and ujt € [wj, tt], j = 3, 4. For example, wi and W3 are the limits for 
contribution of quasiparticle jets along X axis and W2 and W4 those for jets along the Y axis. They can be expressed 
as 



1^1,2= cos [( 1) 6l,2 



1], ^;3,4 = cos-M( " -1)63,4-1], (C5) 



for a pure d-wave ga p, W j — cos~^(— a/6), i = 1, 2, 3, 4. The previously given expressions for Mj{'ip), j = 2, 3, can be 
then used (see Eq. ( ^.1C| )) to obtain M|, M|, the 7r/2 angular Fourier components of the 2a; and 3a; signals of the 
transverse magnetic moment. 
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FIG. 1. Coordinates and definitions used in the paper. The FS and the energy gap are shown schematically. The crystallo- 
graphic directions a and 6 are indicated. The orthogonal X and Y axes are along the nodal directions of the pure d-wave gap. 



The d + s nodal directions, labeled 1,2,3,4 are shifted by an angle ±i/ (see Eq. (2.9)) from their As = values. The applied 
magnetic field, H^, forms an angle ip with the +y direction. 



FIG. 2. Harmonics, Mj, j — 2,3, at 2ui, 3ui of the normalized transverse magnetic moment angular amplitude (Eq. (2. IS)), 



plotted as a function of the field ratio, R = hdc/hac- Normalization is taken so that in the static case, cu ~> 0, Mo is unity. 



FIG. 3. Harmonics, Tj, j = 3,4, at Sui, Auo of the normalized magnetic torque amplitude, (Eq. (2.2C)), as a function of 
R = hdc/hac- Results are shown for gaps with lines of nodes. Normalization is as in Fig. ^ 



FIG. 4. Angular dependence of M3{tf}), (see Eq. 3.9), the 3uj harmonic of the normalized transverse magnetic moment, A4±, 
for various admixtures of s-wave in an d + is energy gap. Results are shown at ^ = hac = 0.05 and normalized (as in Fig. 
^ so that the maximum oi A4± for a pure d-wave gap is unity. Lines are labeled by the ratio of s- and d-wave amplitudes, 
S = As/ Ad. 



FIG. 5. Time dependence in a d-f is state of A^±(V',i) (Eq. (3.11)) and its contributions from jets along the X {A4±{'il),t)x) 
and Y {M±{ip,t)Y) axes. Normalization and field are as in Fig. El Results shown are at fixed angle Tpo — 25" for which 
(Fig. 1^ the signal for S — 0.02 is substantially enhanced as compared to the pure d-wave case. Curves are labeled by the 
axes describing the contribution from a particular jet, and the value of 5. Also shown is cosSujt, which enters the definition of 
M^ixp) in Eq. (3.9) and accounts for part of the enhancement. 

FIG. 6. Origin of the enhancement of Mz{Tp) (the Slj harmonic oi M±_) at h = hac = 0.05, and ip = ipo = 25" (recall Fig. ^) 
with the admixture of s-wave component in the d + is gap. The curves, which represent the product of the normalized magnetic 
moment A4± and cos{3ujt), (see Fig. |^ are labeled by the corresponding value of S. 

FIG. 7. Field ratio, R = hdc/hac, depe nden ce of (the angular 7r/2 Fourier component of the 2lj signal of the normalized 
transverse magnetic moment M± , Eq. (pTCl)), for various admixtures of s-wave in a d 4- is energy gap. Results are given at 
h = 0.05. The normalization is as in Fig.^ Curves are labeled by the ratio S — As/Ad with the same values as in Fig. ^ 

FIG. 8. Field ratio, R, dependence of the angular 7r/2 Fourier component, M3 , of the 3uj signal of the normalized transverse 
magnetic moment, M±. Results shown are for various admixtures of s-wave in a d -I- is state at h = 0.05. Normalization and 
labeling are the same as in Fig. ^. 

FIG. 9. Dependence of the angular it/2 Fourier component, M|, of the 2u! nonlinear magnetic moment signal M± (recall 
Fig. ^) on the admixtures of s-wave component in a d + is gap. Results are shown at h — 0.05 and normalized as in Fig. ^. 
Curves plotted are labeled by the corresponding value of R. 

FIG. 10. Dependence of the angular 7r/2 Fourier component of the 3uj signal, M3 , of A4± (compare with Fig. ^) on the 
admixture of s-wave component in a d -f is state. Results are shown at h = 0.05 and normalized as in Fig. ^. Curves plotted 
are labeled by their value of R. 
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Fig. 1: I. Zutic and O.T. Vails 



Fig. 2: I. Zutic and O. T. Vails, Low Frequency 




Fig. 3: I. Zutic and O.T. Vails, Low frequency 
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Fig. 4: I Zutic and O.T. Vails, Low Frequency 




Fig. 5 : I. Zutic and O. T. Vails, Low Frequency 




Fig. 7: I. Zutic and O.T. Vails, Low frequency 




Fig. 8: I. Zutic and O.T. Vails, Low Frequency 




Fig. 9: I. Zutic and O. T. Vails, Low Frequency 




